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ABSTRACT 

Hydrodynamic loads induced on large underwater structures 
by impulsive motion may be significant design factors. Such 
loads may be induced by earthquake excitation or may result 
from acceleration produced while lifting an object in the 
sea. 

The theoretical approach to the calculation of these 
loads is outlined and numerical results are presented for 
several submerged configurations including a vertical cone, 
a sphere, and a vertical circular cylinder. 

Numerical results for these submerged structures are 
presented in the form of a dimensionless load parameter or 
added mass coefficient. Results corresponding to a number 
of different water depths are presented to show the rather 
sizable effect of the relative water depth on the hydro- 
dynamic force. It is shown that the effect of the free 
water surface is to reduce the hydrodynamic loads in 
comparison to the corresponding infinite depth values. 

Experimental results obtained by vibration testing are 
presented for a submerged sphere, a cone and a vertical 
cylinder. These results show excellent agreement with the 
theoretical results. 
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NOMENCLATURE 



Symbol 

a 



d 

f(5 1 »5 2 »5 3 ) 

& 



g 

G 

h 



£ 

ft 

->• 

n 

P 

P 



S 

t 

3 



^1 * X 2 
u 



Description 

= characteristic dimension of the body 
= added mass coefficient 
= depth of submergence 
= source strength function 
= force vector 
= acceleration of gravity 
= Green's Function 
= depth of water 

= modified Bessel function of the second kind 
of order n 

.= distance of line of action above bottom 
= moment vector 
= unit normal vector 
= dynamic pressure, dimensionless 
= dynamic pressure 
= fluid velocity vector 
= surface of body 
= time 

= velocity vector of object 
= coordinates, dimensionless 
=o 2 a/g 



£-^>£2*^3 



= coordinates of a point on the surface of 
the object, dimensionless 



p 



= fluid density 



Symbol 

a 

4 > 

* 

3 



Description 
= frequency 

= velocity potential, dimensionless 
= velocity potential 
= angular vector velocity of body 
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I. INTRODUCTION 



When an object is accelerated through a fluid there 
are, in general, two types of forces that are recognized, 
one being a drag component and a second, the inertial 
component. The nature of the flow produced by the time 
dependent motion of an immersed rigid object and the relative 
contribution of these two components of force is generally 
considered to be strongly dependent on the amplitude of the 
relative fluid motion in comparison to the characteristic 
lineal dimension of the object. For example, Keulegan and 
Carpenter [Ref. 1] found that, for harmonic motion of the 
fluid past a. fixed circular cylinder, major flow separation 
did not occur and the forces were well represented by 
potential flow values provided the amplitude of the motion 
was less than about a half diameter. If the amplitude 
is small enough the fluid does not move in one direction 
far enough for separation and wake development to occur 
and, accordingly, an inviscid flow analysis is a valid 
mathematical model. 

Accelerated motion of a solid body in a fluid occurs 
for example when a bottom mounted structure is excited by 
an earthquake. Other applications include the accelera- 
tion of ships and accelerated motions of objects being 
lifted in the ocean environment. Southerland [Ref. 2] has 
described the effect of added mass on the mechanics of 
lifting such loads as a submarine rescue chamber. 
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John [Ref. 3] has shown that the velocity potential 
associated with the harmonic motion of a rigid object 
immersed in a fluid with a free surface may be represented 
by a surface distribution of point wave sources for arbi- 
trary values of frequency. The source distribution (or 
weighting) function is obtained in this formulation from 
the solution of an integral equation resulting from the 
kinematic boundary condition applied on the immersed surface. 
Using this technique, numerical results have been obtained, 
for example, for a semiellipsoid oscillating in a free sur- 
face in water of infinite depth by Kim [Ref. 4], the 
intended application being to surface ships. Numerical 
results for a bottom mounted hemisphere have been presented 
by Garrison and Seetharama Rao [Ref. 5 ] and results for 
several axisymmetric configurations by Milgram and Halkyard 
[Ref. 6], Garrison and Chow [Ref. 7 ] have recently presented 
both theoretical and experimental results for wave forces 
acting on a fixed submerged oil storage tank. 

In the work cited, the hydrodynamic force coefficients 
are represented as a function of a dimensionless frequency 

parameter (or equivalently, a dimensionless wave length 

2 — 

parameter) of the form v = a a where a denotes the 

_ S 

frequency of the motion, a denotes the characteristic lineal 
dimension of the submerged object and g denotes the gravi- 
tational constant. This dimensionless parameter is essentially 

a form of the Froude number. In all cases, however, the 
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theory has a fundamental limitation. The source potential 
oscillates with a wave length inversely proportional to 
the dimensionless frequency of the motion v, and consequently, 
a reasonable partition size of the immersed surface limits 
the validity of the numerical scheme to small to moderate 
values of frequency. Specifically, it has been found 
using the Haskind's relations [Ref. 8] and an energy balance 
for purposes of checking accuracy, that results can be 
obtained with reasonable accuracy up to about v = 2.0. 

Beyond this, the accuracy of the numerical results becomes 
questionable . 

Although the range v = 0 to about 2.0 includes most 
applications in ship hydrodynamics and wave forces acting 
on large structures, the frequencies encountered, for 
example, in earthquake excitation corresponds to rather 
large values of v in the case of large structures. Accord- 
ingly, the development of an alternate method which includes 
both the effect of the free surface and bottom on the 
hydrodynamic forces and is valid for large values of the 
dimensionless frequency v is the primary objective of the 
present thesis. Therefore an attempt has been made herein 
to obtain numerical results for several different objects 
for the asymptotic case of very large values of v, i.e., 
fpr the infinite Proude number case. 

High frequency added mass coefficients for two-dimensional 
shapes oscillating on a free surface have been presented 
in a series of papers by Landweber et al . [Ref. 9, 10, 1, 
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and 12]. The motivation for this work was to provide 
hydrodynamic coefficients for two-dimensional ship forms. 

For three-dimensional objects the only results known 
are those of Waugh and Ellis [Ref. 13] for the case of a 
sphere in infinite depth water and Jacobsen [Ref. 14] 
for the case of the cylindrical pier in various water 
depths. In the former case, the theoretical method consisted 
of using successive image doublets to account for the effect 
of ’ the free surface. However, the case of vertical motion 
only was studied and their method was limited to a sphere. 

In the latter case, the solution was represented by a series 
of Bessel functions and is applicable to a vertical circular 
cylinder, extending between the bottom and free surface. 
However, the numerical computations were made using tables 
of Bessel functions with complex arguments which have 
proven misleading. In fact, Jacobsen’s numerical results 
are in error in several instances and therefore, his 
analysis as well as the corrected formulae are presented 
herein. 

In this thesis the theoretical formulation of the 
impulsive hydrodynamic forces acting on a rigid object of 
arbitrary shape is presented. The effect of the bottom 
and free surface is taken into consideration and the force 
coefficient is shown to depend on a frequency of oscillation 
parameter. The problem is first formulated for arbitrary 
values of the frequency and the asymptotic form of the 
solution for both large and small values of the frequer .y 
is discussed. A computer method for numerical evaluation 
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of the force coefficient for submerged objects of arbitrary 
shape is discussed and numerical results are presented 
for a number of specific geometric shapes including a 
circular cylinder, cone, and sphere. High frequency 
experimental results are presented for a vertical circular 
cylinder, sphere, and cone. These experimental results 
are compared with the theoretical results. 
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II. THEORETICAL CONSIDERATIONS 



A. BOUNDARY VALUE PROBLEM 

The geometry involved in the problem under consideration 
is shown in Figure 1. A rigid object of arbitrary shape 
having characteristic lineal dimension a is submerged to a 
depth d in water of depth h and caused to move in a small 
amplitude harmonic motion. The instantaneous motion of the 
rigid object is specified by the lineal velocity 5 and 
angular velocity Q . The object may intersect the bottom 
or free surface. 

Assuming an unseparated, incompressible and irrotational 

->■ 

flow, a velocity potential <f(x,t) defined by 

q = V4> (1) 

where q denotes the fluid velocity vector, may be introduced 
and must satisfy the continuity equation, 

V 2 $ = 0 (2) 

within the fluid region. The bars over the spatial variables 
denote dimensional quantities and x = ix^ + jx 2 + Sx^. 

On the rigid bottom surface described by x 2 = -h the 
fluid velocity in the vertical direction must vanish and, 
accordingly, the corresponding kinematic boundary condition 
on the velocity potential is 



9 $ 





,-h,x 3 ,t) 



0 



( 3 ) 
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The velocity of a point on the immersed surface is given 
by 

^ = ti + q x r (4 ) 

where r denotes the position vector as shown in Fig. 1. The 
appropriate kinematic and dynamic boundary condition on the 
immersed surface, therefore, takes the form 



lit*’*) = t-S (5) 

3n 

where n = tn^ + j"n 2 + icn^ denotes the unit normal vector on 
the surface of the object and is directed outward into the 
fluid region. 

Assuming the amplitude of the motion of the object to 
be small, the kinematic and dynamic boundary condition applied 
at the free surface may be linearized to give the well-known 
free surface boundary condition: 

3 $ (x^,0,x^,t ) 1_ 3 $ ( Xj. , 0 , , t ) _ q (6) 

3X 2 S 3t 2 

Furthermore, if the object oscillates with frequency a, the 
time dependence of the velocity potential may be separated 
and written as Re[<}> (x^ ,x 2 ,x^ )e”^ at ] where Re denotes the real 
part. In which case Eq. (6) becomes 

|i_(x 1 ,0,x 3 ) _ of a 4 ( Xij o,x 3 ) =0 (7) 

where the unbarred quantities represent the coordinates * ade 
dimensionless with a. 
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Fig. 1 - Definition of Geometry 
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In general it is apparent from Eq. (7) that the solution 

to the boundary value problem established must depend on the 
2 — 

dimensionless frequency parameter a a/g. However, there are 

two extreme cases where the results become independent of the 

2 _ 

frequency. These correspond to the case where a a/g 00 in 
which case Eq. (7) becomes 

♦ (x^O.x^) = 0 • (8) 

2— 

For the other extreme where a a/g 0 the free surface appears 
as a rigid boundary as Eq. (7) becomes 



3 <|> ' 
3x, 



( X-i jO ,Xo) 



3=0 



(9) 



If consideration is restricted to either of the boundary 
conditions Eq. (8) or Eq. (9) the velocity potential may be 
written as the sum 

<f>(x,t) = W^t) * 1 (x), i = 1,2, ...6 (10) 



where each term in the sum represents the potential associated 
with the i component of body motion. (i = 1,2,3 refers to 
linear motion components in the x^,x 2 ,x^ directions, 
respectively, and i = 4,5 3 6 refers to the angular motion 
components about the x-^,X 2 ,X 2 axes, respectively.) The time 
dependent function W is defined as 



W ± (t) = U ± (t) a 



^i+3 ( t ) = a 




i = 1,2,3 



( 11 ) 
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in which the symbols and denote the components of the 
lineal and angular velocity vectors of the body, 0 and 
respectively. The potentials <j>^ are dimensionless quantities. 

It is important to note that the application of either of 
the two boundary conditions, Eqs. (8) and (9), makes the 
problem homogeneous in time. Consequently, the assumption 
of harmonic motion may be relaxed and the functions W^(t) may 
be allowed to be arbitrary functions. 

Substituting Eq. (10) into Eqs. (2), (3), (5) and (8) or 
(9) and introducing the dimensionless variables defined as 

x^ = x^/a, h = h/a, d = d/a (12) 

the boundary value problem for (i = 1,2, ...6) may be 
specified as follows: 

2 

7 4> i (x 1 ,x 2 ,x^) = 0 

< M x 1 ,0,x 3^ = °» or = o, (c-*- 0 ) 



where 



9(^ 1 (x 1 ,-h,x 3 ) 



3* 1 (x 1 ,x 2 ,x 3 ) 

9n 



0 (13 a-d) 



m^ (on the immersed surface) 



hi > i 1,2,3 
m., = (x 0 +d)n 0 - x 0 n 0 

42332 (1H a-d) 

m,- = x 0 n n - x,n_ 
b 3 1 13 

mg = x 1 n 2 - (x 2 + d)n 1 



19 



The boundary conditions (13b) represent two different 
possibilities on the free surface, the first being the high 
frequency condition and the second the rigid boundary (or zero 
frequency) condition. In as much as the method of solution 
involving either boundary condition is similar, both solutions 
are developed simultaneously in the following. 

B. GREEN’S FUNCTION SOLUTION 

The solution to the boundary value problem (13) may be 
carried out by use of a Green’s function. The Green's 
function represents the potential for a point source and 
these sources are distributed over the immersed surface 
according to the distribution function f(t). According to 
this concept the potential <f> may be written as 



where (?) denotes a point on the immersed surface, G denotes 



element made dimensionless with the characteristic body 
dimension a. 

In order that Eq. (15) represent a solution to Eqs. (13) 
it is necessary that G satisfy the equation 




(15) 



— —2 

the Green’s function and dS = dS/a denotes the surface area 



v G(x;5) = 6(x-5) 



( 16 ) 



where 6(x) denotes the Dirac delta function, as well as 
boundary conditions (13c) and (13b). 
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The Green's function satisfying Eq. ( 16 ) as well as 
Eq. (13b) and (13c) is obtained by use of successive images 
of sources as depicted in Fig. 2 and is given by 



G(x;t) 



1 

R 



[(-!>" 

1 n=l 





) + (-D 



n+1 




+ 




(17) 



The corresponding form of G satisfying the rigid boundary 
condition Eq. (13c) (a = 0 on x 2 - 0) is also obtained by 
successive images and is similar to Eq. (17) except that all 
of the sources have positive strengths. For this case, the 
Green's function takes the form 



R R. 



0 



R 



Rn 



R, 



R. 



n 



n 



(R 2 

n 


+ + 
R|i 

n n 


R 5 ’ 
J n 


and 


( 18 ) 




/(x i 


- + (x 3 


-* 3 > 2 


m 


^ x 2 ” ? 2^ 




[J~+ 


( x 2 + 5 2 ) 2 




'j ♦ 


(x 2 - 2nh - 


?2 > 2 



( 18 ) 



(x 2 + 2nh + 



(19) 
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22 



Rn -~\J + (x 2 + 2nh - Cp) 2 

n 

R 5 = A/w 2 + (x 2 - 2nh + ? 2 ) 2 

Since G satisfies Eq. (13a-d) the potential as given by 
Eq. (15) must also satisfy these conditions. Thus, it remains 
to select the function f such that Eq. (13d) is satisfied. 
Application of this boundary condition yields the following 
Fredholm integral equation of the second kind 

-f.(x) + ± f f f.(|) dS = 2m 1 (x) , i = 1,2, ...6 (20) 

s 

which is to be solved for f(S). This integral equation is to 
be satisfied at all points x on the immersed surface. 

Because of the form of G, Eq. (20) is rather complex and 
it is therefore necessary to carry out the solution to the 
integral equation numerically. Proceeding in this direction 
the immersed surface may be partitioned into N area elements 
of size AS^ and the integral in Eq. (20) written as the sum 

-f + f a -2m , n = 1,2, ...6 (21) 

n i n. ij n jL * * * 

where i,j = 1,2,...N and n refers to the six possible modes 
of motion of the object. The transformation from the integral 
equation to a summation is possible since f is a well-behaved 
function. It follows from Eq. (20) that the coefficient matrix 
is defined by 



23 



dS 



( 22 ) 



“ij 



ijf |G(x i; t) 

2tt J -£ s 9n 



The function 9G/9n in Eq. (22) is obtained in a straight- 
forward manner by differentiation of Eq. (17) or (18). 

Once the coefficient matrix a^. is obtained from Eq. (22) 
a standard digital computer subroutine may be used to invert 
the matrix equation, Eq. (21), to obtain the source strength 
f at points on the immersed surface. For purposes of evalua- 
tion of G and 9G/9n by use of (17) and (18) it was found that 
about 15 terms were adequate. The data presented herein, 
however, were obtained using 30 terms. 



C. HYDRODYNAMIC FORCES AND MOMENTS 

The dynamic fluid pressure is obtained from the linearized 
form of Bernoulli's equation as 



P = 




(23) 



where the velocity squared term has been neglected. 

t h 

Corresponding to the l component of body motion, the pressure 
is obtained by use of Eq. (10) and (23) as 

= -pw^<j>^ (no sum) i = 1,2,... 6 (24) 

The pressure may further be written in dimensionless coeffi- 
cient form as 

+ M*) 

P 1 (x) = = — <j> i ( x) , i = 1,2, ...6 (25) 

pw A 



i 
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The hydrodynamic force acting on the object is obtained 
by carrying out the integration of the pressure over the 
immersed surface. Using Eq. (24) the force vector associated 
with the j component of body motion is 

» -a 2 JJ* P^.ndS, j = 1,2,... 6 (26) 

The corresponding result for the moment vector is 

P j r x n dS, j = 1,2,... 6 (27) 

Using Eq. (24) as well as the definition of m^ as given in 

t h 

Eq. (14), the i component of force associated with the 
1 1*1 

j component of acceleration is given in dimensionless 
coefficient form as 



Fii n P 

A i1 = “=2 Vi dS ’ 

j Pa 2 w JJ S J 



i = 1,2,3 
j = 1,2,. . .6 



( 28 ) 



and the corresponding expression for moment is 



A 



ij 




-3. 

paw. 




<J> .m. dS , 
J i » 



i = 4,5,6 
j = 1,2,. . .6 



(29) 



where w is defined by Eq. (11). 

Equations (28) and (29) define the inertia (or added mass) 
tensor having, in general, 36 elements. The dimensionless 
coefficients A^. defined by Eq. (28) represent the dimension- 
less 'force coefficient and Eq. (29) represents the dimension- 
less moment coefficient. These coefficients are generally 
referred to as added mass (or moment of inertia) coeff .ents. 
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Although these coefficients are made dimensionless by use of 
a in the definitions (28) and (29), there are cases, depending 
on the body geometry, where other length scales are more 
appropriate. It is also standard practice in many cases to 
use the displaced fluid volume. However, in all numerical 
results presented herein the definition of is specified 

on the figure so that although the definitions vary, no 
confusion should result. 

For purposes of numerical evaluation of the coefficients 
given in Eqs. (28) and (29) the integrals may be written as 
the summation, where the integral sign has been replaced by 
the sum on the k index. 

A, = (<J>,m ) AS, , 1,J * = 1 ’ 2 ’*** 6 ( 3 o) 

0 1 K K k = 1,2, . . .N 

However, in order to apply Eq. (30) it is necessary to first 
evaluate <f>^. at all of the nodal points on the immersed 
surface. For this purpose the surface integration in Eq. (15) 
is also converted to the summation 






j = 1 , 2, ...6 

k,«, = 1,2, .. .N 



(3D 



where the coefficient matrix is defined as 



k£ 




G(x k ;?) 



dS 



( 32 ) 



Once the coefficient matrix is evaluated by Eq. (32), 
the summation indicated in Eq. (31) may be carried out f 



26 



determine <f>^ at all nodal points on the immersed surface. 

The coefficients are then obtained by use of Eq. (30). 

Although there are 36 possible added mass coefficients 

defined by Eqs. ( 28 ) and (29), many are zero on account of 

symmetry and most off-diagonal elements are small. Also, 

even for a body of arbitrary shape A^. is symmetrical so that 

at most only 21 of the coefficients need be evaluated. That 
d 

is, Eq. (13e) along with Eq. (28) and (29) gives 

A iJ 'Us +3 m i ds =JJ S *J ^ ds (33 > 



By use of the Green's theorem applied to the fluid region 
external to the body, it can be shown that 




(34) 



since on the free surface either ♦ = 0 or 9<}>/9n = 0 (depending 
on the boundary condition being enforced) and 9<(>/9n = 0 on 
the bottom. Therefore, in view of Eq. (34) and (33) it is 
apparent that the added mass tensor is symmetrical. 




(35) 



D. JACOBSEN'S SOLUTION 

At this juncture it is appropriate to comment on the 
analytical solution presented by Jacobsen [Ref. 2] for the 
special case of a vertical circular 'cylinder extending between 
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a rigid bottom and a free surface. This solution represents 
the only known analytical result for the problem under 
consideration and, therefore, is of considerable value as a 
point of comparison for the numerical results presented 
herein. However, it was found that a number of errors were 
present in Jacobsen’s work and it was necessary to re-develop 
and re-evaluate the equations. The expression for the added 
mass coefficient associated with horizontal acceleration in 
terms of the notation of the present paper and in terms of 
Bessel functions having real as opposed to imaginary argu- 
ments appears as 



A 



11 



l6h \ 

_3 Z. , 

n=l,3,5 







(36) 



where denotes the added mass coefficient made dimension- 

less with the displaced volume, the factor B n is defined as 



B = 
n 



l/n_ 

V§K> + K7 k 1 ( 5K> 



(37) 



and K and K.. denote modified Bessel functions of the second 
o 1 

kind of order zero and one, respectively. This expression 
agrees with Jacobsen's corresponding expression with the 
exception of a factor of 2.0 and the type of Bessel functions 
involved. It was found that more serious sign errors occurred 
in Jacobsen’s expression for the line of action of the hori- 
zontal hydrodynamic force and, accordingly, his numerical 
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results are also in error. The correct expression for the 
elevation of the line of action (or center of added mass) in 
terms of the notation of the present paper is 




A = n=l,3,5 



2 



(38) 




where £ = T/h, T being the elevation of the line of action 
above the bottom and h the water depth. 

For completeness the expression for the dimensionless 
pressure on the surface of the cylinder is given in terms of 
the notation of the present paper as 



where 9 denotes the angle measured with respect to the x-axis 
in the x-z plane. 

Equations (36) and (37) were programmed for numerical 
evaluation on a digital computer so that an adequate number 
of terms could be evaluated without concern over numerical 
error. These results are plotted in Figs. 8 and 9 for 
comparison with the source distribution results developed 



P 




(39) 
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herein. It was noted that numerical results obtained from 
figures presented by Jacobsen [Ref. 2] are in error when 
compared to numerical results obtained from Eqs. (36) and 
(38). 



« 
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III. EXPERIMENTAL APPROACH 



Experiments were conducted to determine the added mass 
coefficient of a circular cylinder, a cone and a sphere in 
proximity of a free surface and rigid bottom. The basic 
experimental approach consisted of mounting the test 
object on flexible beams and measuring its natural frequency 
of vibration in air and water. The reduction in natural 
frequency when placed in water was then related to the 
added mass of the water. 

The test rig was used to determine the added mass 
coefficient corresponding to the horizontal motion is shown 
in Pig. 3 and 4. The sphere was mounted in this test rig 
on a thin 6 inch long by 2 inch chord vertical strut. Two 
large steel channels were placed perpendicular to each 
other and welded to a 1/2-inch steel plate to form a rigid 
support for the spring-mass system. These channels were 
supported at four points on the top edge of a shallow, water- 
tight box (24" x 70" x 48"). The cylinder, cone, and sphere, 
mounted on a pair of flexible beams, were bolted to this 
supporting structure as shown in Pig. 3. (The cylinder 
is shown. ) 

Four strain gages were mounted on the upper end of 
one of the flexible beams and connected into a Wheatstone 
bridge. The bridge was then connected to a carrier pre- 
amplifier and the frequency of the bridge output was measured 
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Steel channel 




Fig. 3 - Experimental equipment. 
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Pig. 4 - Experimental apparatus. 



by use of a Monsanto programmable counter-timer and in some 
cases by the use of a Visicorder. A ten second duration 
was set on the counter-timer and, depending on the depth of 
water and test object, the readings ranged from about 110 
to 180 cycles. The time duration programmed on the counter 
was highly accurate so that the major error in the process 
was considered to result from reading the number of cycles 
occurring in the 10 second interval to integer values. 
However, since the starting of the interval was random, 
this error was minimized by repeating each test run six 
times and recording the average reading. It was noted that 
the spread in readings from the counter was never greater 
than one cycle in the ten second interval for a given 
configuration. 

A _ calibration curve for the system was developed by 
plucking the mass (test object) in the air with the 
addition of pre-measured weights attached to the test 
object. As a result it was unnecessary to rely on the 
theoretical natural frequency equation for the spring- 
mass system. In general, the test procedure was considered 
to be highly accurate and good repeatability was obtained 
in all cases. The calibration curves are shown in Appendix 
C. 

The same procedure was used for the determination of the 
added mass coefficients for the sphere in vertical motion 
but a different apparatus and tank was employed. The 
parallel beam mount which carried the rod with attached 
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inch diameter plexiglass sphere is shown in Fig. 
6. The tank used in this experiment was 4.0 ft. 
represented, for practical purposes, an infinite 
case. The sphere was oscillated in the vertical 



5 and 
deep and 
depth 
direction. 
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Fig. 5 - Experimental apparatus used in heave tests. 
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IV. DISCUSSION OF RESULTS 



A. CIRCULAR CYLINDER 

Theoretical and experimental values of added mass 
coefficients have been generated for a number of configurations 
and are plotted in Figs. 7-20 and tabulated in Tables 1- 
18. These include results for a vertical circular cylinder, 
a sphere, a hemispherical bottom-mounted vessel, and a 
conically-shaped object. 

The theoretical and experimental results for the added 
mass coefficient for a vertical circular cylinder are shown 
in Figs. 7-12 and Tables 1-7 • The calibration curves are 
presented in Fig. 7* Two different cylinders of the same 
overall dimensions and materials were used. They produced 
calibration curves which were slightly different -result s 
since their natural frequencies of vibration were different 
and their positions in the test rig could not be duplicated 
exactly . 

The added mass coefficient corresponding to the case 
of the horizontal acceleration as based on the cylinder 
radius cubed is presented in Fig. 8, both with a free sur- 
face and with a rigid upper boundary. Results for two 
different grid sizes are shown, 120 and 192 effective nodal 
points, and little difference was found to exist. It can 
be seen that the computer results for the "rigid boundary" 
case closely follow the straight line relationship (uh) 
representing the classical two-dimensional result. The 
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results corresponding to the "free surface" case show that 
the effect of the free surface is to reduce the added mass 
coefficient and for larger values of h the results appear 
to be closely approximated by the straight line relation- 
ship = (irh - 3-0). 

The same type of result in a slightly different form is 
presented in Fig. 9- In this figure the added mass coeffi- 
cient made dimensionless in the more common manner, by use 

— 2 — 

of the displace volume Tra h , is plotted for two different 
grid sizes along with Jacobsen's results, Eq. (36). Also, 
the experimental results obtained from vibration testing 
are presented for comparison. 

The results show that the effect of the free surface 
is to reduce the added mass coefficients and agreement 
between the experimental and theoretical results is- shown 
to be excellent. On the same figure the computer results 
for the rigid boundary case are presented. The classical 
result for two-dimensional flow gives A^ = 1.0 and, 
therefore, the deviation from this value gives the percentage 
numerical error in the computer results directly. The 
figure indicates a maximum of about 3-0 per cent error at 
the extreme low values of h. 

There is no doubt that accuracy of the computer results 
could be improved by means other than the most obvious: 
increasing the number of area elements. However, any gain 
in precision is generally offset by increased computer 
storage requirements, complexity and computer run tine. The 
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ADDED MASS COEFFICIENT 




Fig. 8 - Added mass coefficient for a 
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ADDED MASS COEFFICIENT 



Rigid Boundary 



0.8 - 



0.6 - 
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Fig. 9 - Added mass coefficient for a vertical 
circular cylinder. 
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present method is considered to be quite simple and versa- 
tile with respect to configuration and yet possesses adequate 
accuracy for most purposes. 

For the case of horizontal motion, the line of action 
of the horizontal force is also a parameter of interest. 

When the upper boundary acts as a "wall", the force acts 
at the center, i.e., at l = .5. However, for the free 
surface case the pressures are reduced near the free surface 
and the line of action is somewhat below center. This result 
is presented in Fig. 10 as obtained using 252 nodal points 
and is compared with Jacobsen's corrected result, Eq. (38). 

The figure shows the agreement to be excellent. 

A second degree of freedom is also possible for the 
vertical cylinder. The cylinder may be allowed to rotate 
about some convenient axis with angular acceleration. In 
which case an added moment of inertia may be defined. 

Numerical results have been generated, using 192 nodal points, 
for this coefficient for rotation about an axis passing 
through the center of the cylinder at its base. These 
theoretical results are presented in Fig. 11. The results 
show that the free surface has the effect of reducing the 
hydrodynamic reaction when compared to the rigid boundary 
case. This is the expected result since the free surface 
tends to reduce the pressure differential near the top of 
the cylinder. 

Finally, in order to show how the maximum pressure on 
the surface of the cylinder varies with depth. Fig. 12 as 
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Fig. 10 - Location of effective line of action of hydrodynamic 
inertia force for a vertical circular cylinder. 



Ul h = h/a 




Fig. -11 - Added moment of inertia coefficient 
for a vertical circular cylinder. 
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PRESSURE COEFFICIENT 



1.0 



rigid boundary 




Fig. 12 - Pressure amplitude on meridian of a 
vertical cylinder. 



been prepared. In this figure the dimensionless pressure 
coefficient as defined by Eq. (25) is plotted for the case 
of horizontal acceleration. It may be noted that the 
classical solution corresponding to the rigid boundary case 
yields the maximum pressure of p^ = 1.0 over the complete 
length. Pig. 11 shows that the pressure corresponding to 
the free surface case begins at zero and appears to approach 
unity as the distance from the free surface and depth increase. 
This is the expected trend since the effect of the free 
surface diminishes at large distances. Eq. (39) is plotted 
on the same figure for comparison and the agreement appears 
to be quite good. 

B. SPHERE 

Results for a sphere accelerating both horizontally and 
vertically are shown in Figs. 13-17 and tables 8-l4. The 
calibration data is presented in Fig. 13 (vertical) and 
Fig. 14 (horizontal). For the case of vertical acceleration 
in deep water. Fig. 15 shows experimental results obtained 
by vibration testing for comparison with theory and the 
agreement appears to be quite good. It may be noted that 
the effect of the free surface is always to decrease the 
added mass coefficients while the proximity of a rigid 
boundary has the opposite effect. All of the computer results 
for the sphere have been generated using 264 nodal points. 

The well-known classical solution for a sphere gives an 
added mass coefficient of .5 for the case of a fluid of 
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infinite extent. Accordingly, as d becomes large, the 
coefficient plotted in Fig. 12 is expected to approach 
this closed form result but approximately a 4.0 per cent 
deviation is noted. This error is due to numerical inaccur- 
acies and approximations used in the computer program. 

The same kind of results are plotted in Fig. 16 as 
generated by the computer, the only difference being the 
finite depth, h = 4.0. The same general trends as in Fig. 

12 are observed when the sphere is near the free surface; 
the rigid boundary results are greater than .5 while the 
free surface results are less than .5. However, as the depth 
is increased and the sphere approaches the bottom, the 
added mass coefficient increases in either case. 

Theoretical and experimental results for the added 
mass coefficient .in horizontal acceleration corresponding 
to a bottom mounted sphere is shown in Fig. 17. The 
experimental results were obtained by use of the test rig 
shown in Figs. 3 and 4. The sphere was mounted on the lower 
plate by use of a 6.0 inch strut and the sphere was adjusted 
to within .02 inches of the tank floor. It may be noted 
that the results show good agreement with theoretical 
results and have the same general trends as previous results 
except for the break off in A^ for the rigid boundary case 
as h becomes small. This reversal in trend is - caused by 
the unwetting of the sphere occurring at values of h less 
t han 2.0. 
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Pig. 16 - Added mass coefficient for a sphere 
for h = 4.0. 
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Added mass coefficient for a cotton. 

mounted sphere. 
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C. HEMISPHERE 

Figure 18 and Table 15 present the added mass coeffi- 
cient associated with a bottom mounted hemispherical object 
as a function of water depth for both the rigid boundary and 
free surface case. The results appear somewhat similar 
to previous ones with the rigid boundary having the effect 
of increasing the added mass coefficient and the free surface 
having the opposite effect. 

D. CONE 

Conically shaped bottom mounted, surface piercing 
structures have been proposed for use in the Arctic on 
account of the ice-breaking capability of the shape. This 
structure shown in Fig. 19 is a conically shaped tower 
similar to those which have been proposed as drilling rig 
configurations for deployment in the Arctic. The computer 
results were generated using 192 effective nodal points 
distributed over the immersed surface. The tabular data 
is listed is Tables 16-18. The calibration curve is pre- 
sented as Fig. 20. 

The theoretical and experimental added mass coefficient 
(based on actual displaced volume) for the conically shaped 
tower is presented in Fig. 19 as a function of water depth 
and these results are in good agreement. The results for 
the rigid boundary case shows coefficients which are much 
greater than those corresponding to the free surface case. 
Moreover, for the free surface case the line of action f 
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Pig. 18 - Added mass coefficient for a hemispherical vessel. 




Fig. 19 - Added mass coefficient and line of action 
of horizontal force for a cone. 
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the hydrodynamic force was found to actually lie below the 
bottom. This rather unexpected result is caused by the 
fact that the pressures on the lower portion contribute to a 
moment opposite that induced by pressures on the upper part. 
In general, as- the object is accelerated horizontally, 
positive pressures are produced on the front side while 
negative pressures occur on the back side. Therefore, since 
the pressure as well as area which contributes to a negative 
moment is much greater than those tending to produce a 
positive moment, the result is a moment tending to tip 
the top of the cone in the direction of the acceleration. 
This moment is equivalent to applying the horizontal force 
at a distance T below the bottom. 
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V. CONCLUSIONS 



A practical numerical analysis has been developed using 
three dimensional sources and image sources to represent 
the flow produced during impulsive motion of rigid struc- 
tures of arbitrary shape immersed in water in proximity to 
a free surface. The bottom and free surface have been taken 
into account and it is shown that the proximity of the free 
surface has a sizable effect and always tends to reduce the 
added mass coefficients in comparison to the rigid boundary 
case. The experimental results obtained by vibration 
testing for a sphere, cylinder, and a cone were found to 
agree well with the theoretical results and the results for 
the cylinder are in agreement with analytical results. 

The results are directly applicable to such practical 
problems as the determination of earthquake loading of rigid 
submerged and semi-submerged bottom mounted structures as 
well as the calculation of required forces to lift objects 
of arbitrary shape. 
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APPENDIX A 



TABLES OP COMPUTER GENERATED DATA 



TABLE 1 

CYLINDER: COMPUTER GENERATED DATA (Horizontal Motion) 

192 Effective Points 



h 


Aii=Fn/ p Ui7ra 2 h 

Free Surface 
Boundary Cond. 


A 11 =F 11 /pU 1 Tra 2 h 

Solid 

Boundary Cond. 


n - * 31 

h A 11 

Free Surface 
Boundary Cond 


10.0 


0.910 


1.000 


0.472 


9-0 


0.896 


0.999 


0.469 


8.0 


0.879 


0.999 


0.465 


7.0 


0.861 


1.001 


0.461 


6.0 


0.837 • 


1.004 


0.455 


5.0 


0.807 


1.001 


0.449 


M . 0 


0.764 


1.014 


0.422 


3.0 


0.699 


1.020 


0.434 


2.0 


0.593 


1.027 


0.425 


1.5 


0.512 


1.031 


0.419 


1.0 


0.402 


1.027 


0.414 


0.5 


0.254 


1.062 


0.412 
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TABLE 2 



CYLINDER: COMPUTER GENERATED DATA (Horizontal Motion) 

120 Effective Points 



A 11 =F 11 /pU 1 ira 2 h 
Free Surface 



A^^=F^-^/pU^Tia 2 h 

Solid 



z - 
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Free Surface 








Boundary Cond. 


Boundary Cond. 


Boundary 


10.0 


0.935 


1.006 


0.476 


9.0 


0.908 


1.002 


0.473 


8.0 


0.889 


1.000 


0.469 


7.0 


0.869 


1.000 


0.465 


6.0 


0.844 


1.003 


0.459 


5.0 


0.813 


1.007 


0.453 


4.0 • 


0.769 


1.014 


0.446 


3.0 


0.704 


1.022 


0.437 


2.0 


0.597 


1.031 


0.427 


1.5 


0.515 


1.036 


0.421 


1.0 


0.404 


1.045 


0.416 


0.5 


0.255 


1.069 


0.410 
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TABLE 3 



CYLINDER: COMPUTER GENERATED DATA (Horizontal Motion) 
252 Effective Points 



h 


A 11 =F 11 /pU 1 Tia 2 h 

Free Surface 
Boundary Cond. 


A 1 ^=F 11 /pU 1 ira 2 h 

Solid 

Boundary Cond . 


A 31 

a = _ ^ — 
h A 1X 

Free Surface 
Boundary Cond 


10.0 


0.907 


0.999 


0.471 


9.0 


0.893 


0.998 


0.468 


8.0 


0.877 


0.998 


0.464 


7.0 


0.838 


0.999 


0.460 


6.0 


0.835 


1.002 


0.455 


5.0 


0.804 


1.006 


0.449 


4.0 


0.761 


1.011 


0.442 


3.0 


0.697 


1.017 


0.433 


2.0 


0.590 


1.022 


0.424 


1.5 


0.510 


1.026 


0.419 


1.0 


0.400 


1.032 


0.414 


0.5 


0.250 


1.051 


0.410 
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TABLE 6 



h 

10.0 

9.0 

8.0 

7.0 

6.0 

5.0 

4.0 

3.0 

2.0 
1.5 
1.0 



COMPUTER GENERATED DATA (Rotational Motion 
about vertical axis) 192 Effective Points 



CYLINDER: 

A 66 =M 66 /(pJ V a2 ^ 3) 

Free Surface 
Boundary Cond. 

0.266 
0.258 
0.248 
0.237 
0.225 
0.209 
0.190 
0 . 166 
0.132 . 

0.110 

0.083 



A 66 =M 66 /(p «3’ a2h3) 

Solid 

Boundary Cond. 

0.322 

0.320 

0.319 

0.317 

0.314 

0.312 

0.308 

0.302 

0.294 

0.288 

0.282 
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TABLE 7 



CYLINDER: 

P, 


COMPUTER GENERATED DATA (Pressure 
Coefficients, Free Surface Boundary 
Condition) 


1 PjUja 


d = d/a 


h = 


0.929 


4.722 


5.0 


0.925 


4.167 


5.0 


0.916 


3.611 


5-0 


0.900 


3.056 


5.0 


0.875 


2.500 


5.0 


0.835 


1.944 


5-0 


0.769 


1.389 


5.0 


0.649 


0.833 


5.0 


0.394 


0.278 


5.0 



2.833 


0.847 


3-0 


2.500 


0.841 


3.0 


2.167 


0.827 


3.0 


1.833 


0.804 


3-0 


1.500 


0.769 


3.0 


1.167 


0.716 


3.0 


0.833 


0.637 


3.0 


0.500 


0.510 


3.0 


0.167 


0.279 


3.0 



60 



TABLE 7 cont . 



0.743 

0:735 

0.719 

0.693 

0.655 

0.601 

0.522 

0.405 

0.208 



0.522 

0.515 

0.500 

0.477 

0.445 

0.401 

0.340 

0.255 

0.120 



d = d/a 


h = 


,889 


2.0 


,667 


2.0 


,444 


2.0 


,222 


2.0 


,000 


2.0 


778 


2.0 


556 


2.0 


333 


2.0 


ill 


2.0 



944 


1.0 


833 


1.0 


722 


1.0 


611 


1.0 


500 


1.0 


389 


1.0 


278 


1.0 


167 


1.0 


056 


1.0 



1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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SPHERE: COMPUTER GENERATED DATA (Infinite Fluid Depth) 
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TABLE 11 



SPHERE: COMPUTER GENERATED DATA (h = 4) 

d A 11 =P 11 /pU 1 ^rra^ A 11 =F 11 /pU 1 |ira 3 A 22 =F 22^ P ^2^H^ k^F^pU^na^ 

Free Surface Solid Free Surface Solid 



Boundary 


Cond. Boundary 


Cond . Boundary Cond. 


Boundary Cond. 


3.0 


0.663 


0.682 


0.898 


0.904 


2.9 


0.595 


0.613 


0.717 


0.725 


2.7 


0.556 


0.576 


0.620 


0.632 


2.5 


0.538 


0.559 


0.578 


0.595 


2.3 


0.526 


0.551 


0.554 


0.578 


2.0 


0.514 


0.547 


0.529 


0.570 


1.8 


0.507 


0.549 


0.513 


0.573 


1.6 


0.498 


0.555 


0.495 


0.585 


1.4 


0.485 


O.566 


0.470 


0.610 


1.2 


0.465 


0.590 


0.433 


0.665 


1.1 


0.451 


0.613 


0.407 


0.725 


1.0 


0.433 


0.682 


0.375 


0.904 


0.5 


0.190 


0.225 


0.300 


0.880 
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TABLE 12 



SPHERE: COMPUTER GENERATED DATA (Horizontal Motion 

Bottom Mounted) 



h 


A ll=Pll /(p Ui4Tra 3/3) 


A 11 F n / (P U 




Free Surface 
Boundary Cond. 


Solid 

Boundary 


6 . 0 


0.669 


0.674 


5.0 


0.667 


0.676 


4.0 


0.663 


0.682 


3.0 


0.648 


0.701 


2.5 


0.625 


0.732 


2.0 


0.555 


0.901 


1.90 


0.529 


0.960 


1.62 


0.440 


0.860 


1.43 


0.371 


0.752 


1.22 


0.292 


0.609 


1.00 


0.211 


0.451 



4ira 3 /3) 

Cond. 



TABLE 15 

HEMISPHERE: COMPUTER GENERATED DATA (Horizontal Motion) 



h 


Aii=Pii /(p Ui 2 TTa 3/3 


A ll =P ll /(p ^l 27ra3/3) 




Free Surface 
Boundary Cond. 


Solid 

Boundary Cond. 


6.0 


0.516 


0.518 


5.0 


0.515 


0.519- 


4.0 


0.514 


0.520 


3.0 


0.510 


0.526 


2.5 


0.505 


0.533 


1.7 


0.481 


0.567 


1.4 


0.455 • 


0.610 


1.2 


0.421 


0.674 


1.0 


0.366 


0.901 
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TABLE 16 



CONE: COMPUTER GENERATED DATA (Horizontal Motion) 



h 



A ii- p n / ‘ >u i* 



Free Surface 
Boundary Cond. 



A ii =p ii /c “i v ' 

Solid 

Boundary Cond. 



6.0 


0 


4.0 


0 


3.0 


0 


2.0 


0 


1.5 


0 


1.25 


0 



885 


0.901 


869 


0.922 


845 


0.960 


766 


1.109 


664 


1.405 


578 


1.857 



h 


1 


i 




Free Surface 
Boundary Cond. 


Solid 

Boundary 


6.0 


0.336 


0.336 


4.0 


0.338 


0.335 


3.0 


0.339 


0.333 


2.0 


0.348 


0.324 


1.5 


0.370 


0.305 


1.25 


0.394 


0.281 



Cond. 
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APPENDIX B 



EXPERIMENTAL DATA 



TABLE 4 

CYLINDER: CALIBRATION DATA (Horizontal Motion) 

Cylinder (A) 

Frequency (HZ) Added Weight (lbs) 



18.00 

16.15 

15.70 
15.30 
14 . 1*5 

14.15 
13.85 

13.80 

13.25 

12.70 

11.80 
11.10 



0.00 

0.793 

1.022 

1.275 

1.814 

2.043 

2.296 

2.308 

2.790 

3.329 

4.360 

5.388 



Cylinder (B) 



18.95 

16.63 

15.27 

14.86 

14.10 

13.80 

13.48 

.12.47 

12.04 

11.67 



0.00 

0.864 

1.561 

1.866 

2.351 

2.591 

2.880 

3.888 

4.373 

4.890 



67 



TABLE 5 



CYLINDER: 


EXPERIMENTAL DATA (Horizontal Motion) 


h 


Cylinder (A) 

Frequency (HZ) 


. _ Added Wt 

A 11 Displaced 


5.30 


12.10 


0.835 


O 

OO 

• 

-=r 


12.50 


0.821 


4.20 


13.00 


0.804 


3.70 


13.48 


0.775 


3.15 


14.10 


0.722 


2.75 


14.55 


0.682 


2.25 


15.30 


0.618 


2.00 


15.66 


0.585 


1.65 


16.10 


0.548 


1.25 


16.75 


0.441 


0.90 


O 

C\J 

• 

r— 

rH 


0.380 


0.40 • 


17.78 


0.221 


0.00 


18.00 


0.000 




Cylinder (B) 




5.30 


12.43 


0.818 


4.80 


12.80 


0.816 


4.20 


13.40 


0.788 


3.70 


13.92 


0.745 


3.15 


14.60 


0.694 


2.75 


15-07 


0.678 
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TABLE 9 



SPHERE: CALIBRATION DATA 

Frequency (HZ) 

19.8 

19.0 

18.6 

18.3 

17.5 

16.7 

16.4 

16.0 

15.4 
15.1 

15.0 

14.6 

14.4 



(Vertical Motion) 

Added Weight (lbs.) 

0.0 

0.3 

0.5 

0.7 

1.0 

1.3 

1.5 

1.7 

2.0 

2.3 

2.5 

2.7 
3.0 
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TABLE 10 



SPHERE: EXPERIMENTAL 


RESULTS (Vertical Motion, Infinite fluid 
Depth) 


h 


. _ Added Weight 

r 22 Displaced Weight 


6.0 


0.514 


5.0 


0.480 


4.0 


0.514 


3-0 


0.504 


2.0 


0.500 


1.7-5 


0.475 


1.5 


0.463 


1.25 


0.416 


1.0 


0.367 
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TABLE 13 

SPHERE: CALIBRATION DATA (Horizontal Motion, Bottom Mounted) 
Frequency (HZ) Added Weight (lbs.) 



14.65 


0.379 


14.35 


0.592 


14.00 


0.828 


13.95 


0.864 


13.65 


1.077 


13.40 


1.313 


13-00 


1.625 


12.75 


1.873 


12.25 


2.390 


11.85 


2.875 


11.45 


3.397 


10.75 


4.405 
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TABLE 14 



SPHERE: EXPERIMENTAL DATA (Hori 



h 


Frequency 

HZ 


3.72 


12.05 


3.40 


12.00 


3.13 


12.05 


2.90 


12.10 


2.67 


12.15 


2.50 


12.20 


2.33 


12.25 


2.17 


12.35 


2; 00 


12.45 


1.83 


12.70 


1.67 


12.90 


1.50 


13.20 


1.33 


13.55 


1.17 


13.85 


1.00 


14.15 


0.00 


15.32 



ontal Motion, Bottom Mounted) 

A _ added wt. 

±1 ~ displaced wt . 

0.644 

0.661 

0.644 

0.624 

0.612 

0.605 

0.585 

0.563 

0.536 

0.465 

0.421 

0.355 

0.287 

0.230 

0.179 

0.00 
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TABLE 17 



CONE: CALIBRATION DATA (Horizontal Motion) 



Frequency (Hz) 


Added Weight (lbs.) 


17.61 


0.418 


16.89 


0.631 


16.50 


0.867 


15.53 


1.352 


15.11 


1.664 


• 14.62 


1.912 
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TABLE 18 



CONE: EXPERIMENTAL DATA (Horizontal Motion) 



h 


Frequency (Hz) 


. _ Added Wt . 

A 11 Displaced 

Wt . 


0.79 


15.32 


0.281 


0.64 


16.03 


0.214 


0.48 


16.89 


0.142 


0.33 


17.79 


0.095 


0.17 


18.66 


0.065 
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CALIBRATION CURVES 
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Fig. 7 - Calibration Curve (CYLINDER) 



ADDED WEIGHT (LBS.) 




Fig. 13 - Calibration curve • 

(Sphere - Vertical Motion) 
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Fig. 1^4 - Calibration curve (Sphere - Horizontal Motion) 
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Fig. 20 - Calibration curve (Cone - Horizontal Motion) 
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